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Abstract 

A (64, 40, 8) subcode of the (64, 42, 8) third-order Reed- Muller code is proposed to NASA 
for high-speed satellite communications. This code can be either used alone or used as an 
inner-code in a concatenated coding system with the NASA standard (255,223,33) Reed- 
Solomon code as the outer code to achieve high performance with reduced decoding com- 
plexity. This Reed-Muller subcode has a relatively simple and parallel trellis structure and 
consequently can be decoded with a group of identical and relatively simple Viterbi decoders 
in parallel to achieve high-speed decoding. In this report, the complexities of various sec- 
tionalized trellis diagrams are analyzed. Based on this analysis, the trellis diagram with the 
smallest overall complexity will be used for the implementation of a high-speed decoder. 
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1. Introduction 


i ■■ 

j- 


A (64,40,8) subcode of the third-order Reed-Muller (RM) code is proposed to NASA 
for high-speed satellite communications. This RM subcode can be used either alone or as an 
inner code of a concatenated coding system with the NASA standard (255,223,33) Reed- 
Solomon (RS) code as the outer code to achieve high performance with reduced decoding 
complexity. It can also be used as a component code in a 3-level bandwidth efficient coded 
8-PSK modulation system together with the (64, 7, 32) first-order and (64, 57, 4) fourth-order 
RM codes. 

The (64, 40, 8) RM subcode has a relatively simple and parallel trellis structure and 
consequently a group of identical and relatively simple Viterbi decoders can be devised to 
process the decoding in parallel. This not only simplifies the decoding complexity but also 
speeds up the decoding process. As a result, a high-speed decoder can be implemented. For 
the AWGN channel using BPSK transmission, this code achieves 5.3 dB coding gain over 
uncoded BPSK at the bit-error rate (BER) of 10 -6 as shown in Figure 1. If this code is used 
as the inner code in a concatenated coding system with the NASA standard (255,223,33) 
RS code as the outer code, the concatenated system will achieve error-free communication 
for signal-to- noise ratio (SNR) greater than 3 dB as shown in Figure 1. It achieves a 7.6 dB 
coding gain over the uncoded BPSK system at the BER of 10 -6 . The (64,40, 8) RM subcode 
has a rate of 0.626 bits/symbol which is higher than that of the NASA standard rate-1/2 
64-state convolutional inner code. In performance, the (64,40,8) RM code achieves a 0.5 dB 
less coding gain than the NASA standard convolutional inner code. Because the (64,40,8) 
RM subcode has a higher rate and can be implemented at a much higher decoding speed, 
it is an attractive alternative to replace the NASA standard convolutional inner code in a 
concatenated coding system with only a small loss in coding gain. 

In this report, we investigate the complexities of various sectionalized trellis diagrams for 
the (64, 40, 8) RM subcode. Particularity, we investigate the complexities of the 4-section, 8- 
section and 16-section trellis diagrams for the (64,40,8) RM subcode. Based on this study, 
we will choose the trellis diagram which has the smallest overall complexity for decoder 
implementation. 
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2. Trellis Complexity of Binary Linear Block Codes 


Any linear block code can be decoded by applying the Viterbi decoding algorithm 
to its minimal trellis using either Euclidean distance or correlation as the distance metric 
[1]. For BPSK transmission over an AWGN channel, the latter has distinct computational 
and implementation advantages. We define the computational complexity of any trellis as 
the number of addition equivalent operations needed to decode a single real valued received 
vector into the maximum likelihood solution (the codeword which has the highest a posteriori 
probability) by using the Viterbi algorithm based on the trellis. In a paper by Kasami et. 
al. [3], the computational complexity of a minimal L-section trellis of a block code has been 
calculated in terms of: (1) the number of real number comparisons, (2) the number of real 
number multiplications and (3) the number of real number additions. In the present case, 
by choosing to maximize the correlation between the received vector and the codewords we 
avoid the use of real number multiplications while decoding a received vector. Furthermore 
the complexity of a real number comparison is roughly the same as that of a real number 
addition. Hence we are justified in expressing the computational complexity of the various 
trellises in terms of the number of real number addition equivalent operations. 


Definition: 1 An L-section trellis diagram for a binary linear (N,K) block code C with 
N = LM is a directed graph T with LA- 1 levels (henceforth called depthsj labelled 
0, M, 2 M, . . . , ML and vertex sets So, Sa/, • • . , Sml (henceforth called state setsj such that 

1 . So — [po] , Sn — Slm — { Pf } where p Q and pp denote the initial and final states 
respectively. 

2. There are edges (henceforth called branches,) connecting states in adjacent state sets 
SiM, *%+i )M: 0 <i < L with each branch labelled by a binary M-tuple and originating 
at a state in S{m- 

3. There is a directed path from the initial state p 0 to the final state pp with label sequence 
Ui , u 2 , . . . , u i if and only if 

(ill o u 2 0 • • • 0 u L ) 

is a codeword in C where o stands for concatenation of two sequences. 

Let C denote an (iV, d tn j n ) linear block code. Let L , M be positive integers such 
that LM — N. Given a L-section trellis, the set of states at the end of each section 
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{So, 5 m, S 2 m, • • • , Slm), we define a sequence {.So? • • • > $lm} called the state 

dimension profile (SDP) of the trellis and given by S{m — log 2 (|5,M|)- The minimal 5- 
section trellis of a code C has the property that every component of its SDP is less than or 
equal to the corresponding component in the SDP of any other proper 5-section trellis for 
C. The maximum among the TV-f 1 components in the SDP of the minimal N - section trellis 
for C is denoted s max (C) and we will denote the maximum of the components in the SDP 
of the minimal 5-section trellis for C as s maXi L(C). For a binary TV-tuple v = (iq, . . . , tqv), 
let Ph,h'[v] denote the ( h ' — h)- tuple . . . ,Vh>) and let 

Ph,h' [51] — {ph,h' [c] • C £ C) 

Let C/^/ be the linear subcode of C consisting of all codewords whose components are all 
zero except for the ( h 1 — h ) components from the ( h -f- l)-th bit position to the Zi'-th bit 
position. 

The computational complexity of any trellis T is denoted A(T) and is composed of two 
components, one as a result of branch metric computations denoted Ab{T) and the other 
due to the state metric computations and comparisons denoted Al{T) yielding AiT) — 

A b (T) + A l (T). 

In a 5-section minimal trellis for a block code, there may be a set of parallel branches 
between two adjacent states. In such a case, we call the entire set of parallel branches a 
composite branch. Each composite branch in the 5th section 1 < i < 5, is made up of 
2 Pi parallel branches where P, is the dimension of the subcode denoted as shown 

in [2]. In the 5th section of a 5-section trellis for a linear block code 1 < i < 5, the number 
of distinct branch metrics that have to be computed is 2 Di where D{ is the dimension of the 
subcode P(i_i)M,iA/(C) and this number is much less than the total number of branches. A 
is the rank of the submatrix formed by M columns from the ((i — I )M + l)-th to the (iM )- th 
column of the generator matrix of the code and is upper bounded by M . For 1 < i < 5, let 
the number of composite branches merging into any state s € Sim be 2 SiM (it is the same for 
any state in 5 {m)- For 0 < i < 5, let the number of composite branches emanating from 
any state s € SiM be 2 A,M , (it is the same for any state in S{m)- &iM is related to as 

follows: 

<W = — s iM 

Therefore the { S{m : 1 < i < 5} can be computed from the SDP of the code and {A,m : 
0 < i < L}. Based on the theory of 5-section trellises [2, 4] it can be shown that 

A( t -_i)M = dim(C(i_i)A/,jv) — dimfCuv/.iv) — Pi 

These dimensions can be easily determined from the trellis oriented generator matrix of the 
code [4]. 
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For the special case of Reed-Muller codes and L ~ 2 l for a positive integer /, since the 
L- section minimal trellis is symmetric about the midpoint [2], i.e., the last X/2 sections form 
a mirror image of the first L/ 2 sections and we have the following. 

$iM ~ *^( L—i)Mi ^ — l — T/2, ^-Ild \ iM — ^ 

The five parameters S{m, DiM, PiM, Km, $iM discussed above are essential to calculating the 
computational complexity of a minimal trellis. 

Consider the operation of a Viterbi decoder using a T-section minimal trellis of a linear 
block code. In order to compute the state metric of a state s E SiM> 1 < * < T, the decoder 
has to pick the minimum among 2 SiM candidates. For 1 < j < 2 6,M , the y -th candidate 
metric is given by the sum of the metric of the j-th composite branch and the metric of 
the originating state of the j- th composite branch. Thus are 2 6iM additions and 2 S,M — 1 
comparisons are required to compute the state metric of s. Hence we have the number of 
addition equivalent operations as 

L 

A l (T) = (2^ M+1 - 1) (2.1) 

»=i 

In the Fth section of the trellis, the decoder has to compute the 2 D ' distinct simple 
branch metrics at the cost of (M — 1) additions per branch metric. If the branch metric 
computations are done in the most parallel manner with 2 Di adders operating simultaneously 
in parallel and independent of each other, a total of 2 Di (M - 1) additions is required. There 
are other slower methods of computing the same set of branch metrics that use lesser number 
of additions. In addition to the distinct simple branch metrics, the branch metric of each 
composite branch needs to be computed. The branch metric of a composite branch is the 
minimum of the branch metrics of the 2 P ' simple branches that form a composite branch. 
Given the branch metrics of all the simple branches, 2 Pi - 1 comparisons are required to 
compute the metric of a composite branch. Therefore, the number of additions required to 
compute the metric of all composite branches is given by the number of distinct composite 
branches times (2 Pi — 1). Based on the above discussion we have, 

L L 

A B (7) = ^ D ' (M — 1) + ^(Distinct composite branches in section-i) ( 2 Pi — 1) (2.2) 

t'=l i=l 

Based on the above discussion, a study of all the possible L-section trellises for the (32, 16, 8) 
Reed-Muller code yielded Table 1. 

From Table 1, it is evident that the computational complexity as well as the maximum 
state space dimension depend on the number of sections in the trellis of the code. For the 
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above code, although A(T) reaches a minimum of 4, 064 when L = 8, the trellis with L = 4 
and A(7~) = 6, 847 is more suitable for implementation because while the former has only 
2 parallel isomorphic trellises with each subtrellis having s max = 128 states the latter has 8 
parallel isomorphic subtrellises with each subtrellis having s max = 8 states. 


3. Trellises with Parallel Structure 

The objective of this section is to show that we can build a trellis for a linear block 
code say C which is a union of certain desired number of parallel isomorphic subtrellises. 
Although this trellis is not minimal , it has a maximum state space dimension less than or 
equal to <s max (C). The conditions under which such a trellis construction is possible and 
an upper bound on the number of such parallel subtrellises is derived. In some cases the 
minimal trellis itself possesses a parallel structure. The number, of such parallel subtrellises 
(if any) in the minimal trellis is derived. 

Suppose we want to build a trellis for linear block code C as a union of 64 parallel 
isomorphic subtrellises. This can be done by first choosing a subcode CT of dimension 6 
less than the dimension of C. Second, we build the trellis for Ck Finally, the trellis for 
C can be obtained as the union of 64 parallel isomorphic subtrellises with each subtrellis 
being isomorphic to the trellis of Ck The resulting trellis for C may not be minimal. We 
prove necessary and sufficient conditions for the existence of a subcode C' of C such that 
the resulting trellis for C has maximum state space dimension s max (C) while retaining the 
property of parallelism. 

Trellises with parallel structure are desirable for the VLSI implementation of high speed 
decoders for block codes. One can envision a decoder consisting of as many parallel processors 
as there are parallel subtrellises in the total trellis, each processor producing a best candidate 
from its subtrellis after processing it with the Viterbi algorithm. Finally it is a simple matter 
of choosing the best codeword among the choices presented by the processors. 

Parallel Trellises with Constraint on Maximum State Space Di- 
mension 

Given a linear block code, several non-isomorphic trellises are possible. The number of 
such non-isomorphic trellises is related to the placement of non-attacking rooks on an upper 
triangular chess board [5]. However there is a unique proper trellis called the minimal trellis 
such that its state space dimension is less than or equal to the state space dimension of any 
other proper trellis at every position. A systematic method for the construction of the 
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minimal trellis using the trellis oriented generator matrix has been presented in [4]. One 
drawback of the minimal trellis is that often it does not have any parallel structure; or the 
parallel subtrellises are themselves too large. In such a case it is desirable to explore the 
possibility of increasing the parallelism of the trellis even if minimality has to be compro- 
mised. A trellis construction with large number of parallel subtrellises for the special case 
of Reed-Muller codes has been presented in [6]. This so called primitive decomposition has 
the drawback that the state space dimension is very large and usually much larger than 
the maximum for the minimal trellis. For many codes the primitive decomposition is quite 
attractive. 

Let G be the trellis oriented generator matrix of an ( N,K ) binary block code C. Let 
r = (r 1} r 2 , . . . ,r N ) be a typical row of G. Then, we define the span of r denoted span( r) 
to be the smallest interval [i,j], 1 < 2 < j IV which contains all the non-zero elements of 
r. For a row r whose span is [i,j] we also define an active span of r denoted aspan( r) as 
IhJ- 1] if i < j and aspan( r) = 4> if i — j. The trellis oriented matrix has the following 
properties: 1) The leading 1 of every row occurs in an earlier position than the leading 1 
of the row below it. 2) The trailing 1 of every row occurs at a different position from the 
trailing 1 of every other row. 3) Any other trellis oriented matrix for C has the same set 
of row spans although the rows themselves may be different [5]. Let 'T be the minimal 
n-section trellis for C. It can be proved that given the trellis oriented generator matrix of a 
code, the state space dimension at any position / is just equal to the number of rows whose 
active span contain l [5, 4]. For example, consider the following trellis oriented generator 
matrix 


/ 1 

i 

1 

1 

0 

0 

0 

0 

r i 

\ 

0 

i 

0 

1 

1 

0 

1 

0 

T2 


0 

0 

1 

1 

1 

1 

0 

0 

T 3 


V 0 

0 

0 

0 

1 

1 

1 

1 

r 4 

) 


for which aspan(yi) = [1,3], aspan(v- 2 ) = [2,6], aspanfa) = [3,5] and aspan^r^ = [5,7]. 
For each /, 0 < / < 8 counting the number of rows which are active at each / for 0 < l < 8 
yields the state dimension profile (SDP) { 0 , 1, 2, 3 , 2 , 3 , 2 , 1, 0 }. For 0 < / < N, let si{C) 
denote the dimension of the /- th state space of C. Let s max (C) be the maximum among the 
state space dimensions. Define the non-empty set 

/max(C) = {l : a,(C) = w(C)} (3.3) 

Suppose we choose a subcode Cfi of C such that dim(C / ) = dim (C) — 1 and the set of coset 
representatives [C/C'] is generated by the single row r £ G. From the above statement 
about s/(C), it is clear that s/(C') = si( C) - 1 for exactly those / where r is active, i.e., 
I £ aspan(r). For other positions / ^ aspan( r) we have si( C') — s/(C). Hence we have the 
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following proposition. 


Proposition: 1 If there exists a row r in the trellis oriented generator matrix G for the 
code C such that aspan(r) D 7 max (C) then, we can form a subcode C' of C generated by 
G Such that -S max (c ) — •Smax(C) 1- UTid 7 max (C ) — 7 max (C). 

In fact 7 max (C') = 7 max (C) U {/ : s/(C) = s max (C) — 1 ,1$ aspan(r)}. Since G is a trellis 
oriented generator matrix, G' = G — {r} is also trellis oriented. We can apply the above 
proposition again to C' if there exists another row r' 6 G' with aspan( r') D 7 max (C'). This 
yields a subcode C with dimension smaller by one and <s max (C) = s max (C ) — 1. If no such row 
r' exists the proposition cannot be applied and the recursion stops. The above proposition 
can be generalized. 

Let 72( C) be the subset of rows of G, the trellis oriented generator matrix of a code C 
consisting of those rows r given by 

R( C) = {r 6 G : aspan( r) D 7 max ( C)} (3.4) 

Let d = 1 72(C) | where \Q\ denotes the cardinality of any finite set Q. 

Proposition: 2 With 72(C) defined as above and d = |72(C)| ; let \ < d' < d. There exists 
a subcode C' of C such that s max (C / ) = s max (C) — d' and dim(C / ) = dim(C) — d! if and 
only if there exists a subset R' C 72(C) consisting of d' rows of 72(C) such that for every 
l satisfying 5/(C) > 5 max (C / ) there exist at least s/(C) — s m ax(C') rows in 72' whose active 
span contain l. The set of coset representatives [C/C'] is generated by R ' . 

Proof: Suppose R' = {r'l, . . . , r'gi] satisfies the conditions in the hypothesis. Since 72' C 
72(C), 7 max (C) C aspan(r'i) for 1 < i < d 1 . Consider the subcode generated by G — 72'. For 
those l 6 7 max (C), we can determine si(C') by counting the number of rows r E (G — 72') 
that are active at that position /. But this number is exactly less than s Inax (C) by 7'. For 
l $ 7 max (C) and satisfying sfi C) > s max (C'), we are assured by the hypothesis that s/(C) 
will be reduced by at least s<(C) — 3 max (C') thus guaranteeing that s max (C') — ^max (C )-<e. 

To prove the converse, let C' be a subcode of C whose dimension is dim(C) — d! and 
satisfying s max (C') = s max (C) — d! . Without loss of generality we may let C' be generated 
G — 72' for some subset 72' of the trellis oriented generator matrix G of C with |72'| = d! . 
Let T be the minimal trellis corresponding to G. Let T' be the minimal trellis for C'. Let 
TV/ (72') be the number of rows r' in 72' such that / E aspan( r'). Then, at every position 
/, 0 < / < TV, we have 


Sl {7) = sfiC') + 7V/(72') > 5,(C) 


(3.5) 



since s/(C) is the smallest possible state space dimension. Therefore 


> s,(C)-s,( C') (3.6) 

Ni(Bf) > a/(C)-w(C') 

For every / at least s/(C) — s max (C') rows of R! are active. Also, for every l E I m ax(C) we 
have Ni(R') > s max (C) — s max (C') = d! . So all the rows r 7 € R' satisfy aspan(r') O ^max(C). 
Thus i? 7 C R( C). 9 

The utility of the above proposition is that it shows how to choose a subcode C 7 of 
C with 5 max (C / ) = s max (C) — dim ([C/C 7 ]), such that one can build a trellis T for C which 
although not minimal has the following properties. 


1. The maximum state space dimension of T is s max (C). 


dim 


fC/C'l 


2. T is the union of 2 i -1 parallel isomorphic subtrellises T t - with each T,- being 
isomorphic to the minimal trellis for Cb 


3. The smallest such subcode has dimension lower bounded by dim(C) — |i?(C)|. i.e., 
the maximum number of parallel subtrellises one can obtain with the constraint that 
the total state space dimension never exceed s max (C) is upper bounded by 2^^ with 
R(C) as defined above. 


The following proposition follows from the above discussion. 

Proposition: 3 The logarithm to the base 2 of the number of parallel isomorphic subtrellises 
in a minimal L-section trellis for a binary (n, k ) linear block code is given by the number of 
rows in its trellis oriented generator matrix whose span contain the integers {M, 2 M, . . . , (L — 
1 )M} where N — LM . 


As an example consider the extended code BCH( 32,21,6) for which the minimal 4- 
section trellis has SDP {0, 7, 9, 7,0). Any trellis oriented generator matrix for this code has 
the following unique set of row spans. 


span( ri) = [1,8] 

span(r 3 ) = [3, 13] 
span{ r 5 ) = [5, 12] 
spanipof) = [7,21] 


span( r 2 ) = [2, 15] 
span( r 4 ) = [4, 14] 
span( r 3 ) = [6, 18] 
span( rg) = [8, 25] 
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span( r 9 ) = [9,16] span(r i0 ) = [10,23] 
5pan(rn) = [11,19] spcm(r 12 ) = [12, 26] 

span( r 13 ) = [13,20] span(ri 4 ) = [14, 22] 

span(r 15 ) = [15,27] span(rie) = [17, 24] 

span( r 17 ) — [18,31] span(r 18 ) = [19, 29] 

span( rig) = [20,30] s/wm(r 2 o) — [21, 28] 

span(r2i) = [25,32] 


4nax(C) = {16} and it can be verified that |i?(C)| = 9. In an attempt to build a trellis 
consisting of 64 parallel subtrellises while satisfying the upper bound of 9 on the maxi- 
mum state space dimension, we let d! — 6. So <s max (C) — d! = 5 max (C / ) = 3. The set 
{/ : si( C) > s max (C')} = {8,16,24}. However, we find that no subset R! of R( C) exists 
satisfying the conditions in proposition 2. Hence we cannot build a trellis consisting of 64 
parallel subtrellises for this code without violating the constraint on the maximum state space 
dimension. If we choose d' = 5 then we can find a subset R' = {r 8 , r 7 , r 8 , r 12 , r 15 } C R(C) 
that satisfies all the conditions in proposition 2. Hence choosing the subcode C / generated 
by G — R' we obtain a trellis T for C consisting of 32 parallel isomorphic subtrellises. Each 
subtrellis is isomorphic to the minimal trellis for C' which has s max (C / ) = 4. 


4. The Trellis with Least Complexity for the (64, 40, 8) 
Code 

Let C denote the RM (64, 42, 8) code and C a (64, 40) subcode of C. We want to examine 
all possible L-section trellises with an aim to find their suitability for decoding. We realize 
that if the trellis to be used for decoding purposes has a parallel structure and is composed 
of a union of M parallel subtrellises, then M smaller Viterbi decoders can be designed to 
process each subtrellis simultaneously in parallel. Therefore the effective complexity is 
just that of a single subtrellis plus the cost of obtaining the final decision from each of the 
M Viterbi decoders. We determine the value of L which minimizes the effective complexity 
with the constraint that the L-section trellis T have a maximum state complexity not greater 
than s maXv L(C). This constraint is reasonable because although one can essentially increase 
the number of parallel subtrellises without bound the hardware complexity of implementing 
the resulting trellis and the cost of combining outputs from subtrellis decoders both increase 
exponentially in the state-space dimension of the code. 

Thus, we will determine the smallest effective complexity denoted H e //(L), that can 
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be achieved with an L-section trellis satisfying the constraint on maximum state dimension 
among the choices L — 1,2, 4, 8, 16, 32, 64 and pick that L which yields the minimum A e fj(L). 
The uninteresting choices of L = 1,2,32,64 were eliminated from consideration by analysis 
which is not presented here. Herein, we consider consider the more interesting cases of 
L — 4, 8 and 16. 

4 . 1 . £ = 4,M = 16 

| 

Let C 0 — (16,15,2), Ci = (16,11,4) C 2 = (16,5,8) be the corresponding Reed-Muller 
I' codes, G, a generator matrix of C* and G{/j a generator matrix for the set of coset repre- 

sentatives [C t /Cj]. For L — 4 the RM (64,42) code has a minimal trellis (with 16 parallel 
f subtrellises) corresponding to the 2-level squaring construction with a state dimension profile 

| (SDP) {0, 10, 10,10,0} (-s maX)4 (C) = 10) and trellis oriented generator matrix 

G = ( 1 1 1 1 ) ® G 0/1 + 

/I 1 0 0 \ 

0 110 0 Gj/2 t 

V 0 0 1 1 / 

10 0 0 
0 10 0 
0 0 10 
0 0 0 1 

In order to obtain a generator matrix for C, the (64,40) subcode of C, one can delete any 
two of the 64 rows above. The s maX;4 (C) of the resulting code depends on which two rows 
we delete. It is easy to see that in order to have the least s max (C,4) which equals 8 we 
must delete any two of the 4 rows among (1111) 0 G 0 /i obtaining an SDP of {0,8, 8, 8,0} 
(s m ax, 4 (C) = 4). Thus 7 max (C) = {1,2,3} and R{ C) = {(1111) 0 G' 0 /i} where G' 0 /i has 
dimension 2. Since |i?(C)| = 2, we can obtain at most 4 parallel subtrellises in any 4-section 
without exceeding the allowable s maXf4 of 8. Choosing C' as the subcode of C generated by 

/ 1 1 0 0 \ 

G' = 0 110 0 G 4 /2 T 

\ o o i i / 

/ 1 0 0 0 

0 10 0 

0 0 10 

V o o o i 
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we can obtain a trellis for C as the union of 4 parallel isomorphic subtrellises T\, i = 
0, 1, 2, 3. The complexity parameters of the subtrellises is presented in Table 2, The resulting 
computational complexity parameters are 

A B (T'i) = 23,296 and A L (T' { ) = 16,383 (4.9) 

This gives A e //(4) = 23,296 + 16,383 + 3 — 39,682 addition equivalent operations. 

4.2. L = 8, M — 8 

Let Co = (8,8,1), Cl = (8,7,2) C 2 = (8,4,4) C 3 = (8,1,8) be the correspond- 
ing Reed-Muller codes. For L — 8, the RM (64,42,8) code has a minimal trellis (with 
2 parallel subtrellises) corresponding to the 3-level squaring construction with a SDP 
{0, 7, 10, 13, 10, 13, 10, 7, 0} (s maX)S (C) = 13) with trellis oriented generator matrix 

G = ( 1 1 1 1 111 1 rg ) (8) G 0 /i + 

/ 1 1 1 1 0 0 0 0 rj 

0 1 0 1 1 0 1 0 rj 

0 0 1 1 1 1 0 0 r\ 

\ 0 0 0 0 1 1 1 1 r\ 

1 1 0 0 0 0 0 0 

0 1 1 0 0 0 0 0 

0 0 1 1 0 0 0 0 

0 0 0 1 1 0 0 0 

0 0 0 0 1 1 0 0 

0 0 0 0 0 1 1 0 

0 0 0 0 0 0 1 1 

The (64, 40) subcode C with the best SDP is obtained by deleting the rows rj] 0 G 0 /i and 
any one among the three rows r\ 0 Gi/ 2 - This code C has SDP {0,6,8,11,8,11,8,6,0} 
(«s ma x,8(C) = 11) and thus / max (C) = {3,5}. By inspection R(C) = {r\ 0 G!/ 2 ,rJ 0 G' 1/2 } 
yielding |i?(C)| = 5. Thus one can obtain at most 32 parallel subtrellises in any 8-section 
trellis for C without exceeding the maximum allowable state space dimension of s maXi8 (C) = 
11. Let C' denote the (64,35) subcode of C generated by G / = G — R{ C) with SDP 
{0, 6, 6, 6, 3, 6, 6, 6, 0}. 

Hence one can build a trellis for C as the union of 32 parallel isomorphic subtrellises 
T 7 ,-, 0 < i < 32. The complexity parameters of the subtrellises is presented in Table 3. The 
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resulting computational complexity parameters are 

Ab(T'<) = 6,784, and A L (T f { ) = 6, 007 (4.11) 

This gives 4 e //(S) = 6,784 + 6,007 + 31 = 12,822 addition equivalent operations. 

4.3. L = 16, M = 4 

Let C 0 = Ci = (4,4,1), C2 = (4,3,2) C3 = (4,1,4) C 4 = (4,0, 00) be the correspond- 
ing Reed-Muller codes. For L = 16, the RM (64,42,8) code has a trellis oriented generator 
matrix given by 

G = Grm(16,5,8) <8> Gi/2 + Grm(16,11,4) ® C*2/3 + GrM(16,15,2) ® G3/4. (4.12) 

where Grm(ti,m) denotes a trellis oriented generator matrix for the corresponding Reed- 
Muller code. For L = 16, the RM (64,42,8) code has a minimal trellis (with no 
parallel subtrellises) corresponding to the 4-level squaring construction with a SDP 
{0,4,7,10,10,13,13,13,10,13,13,13,10,10,7,4,0} (s max , 16 (C) = 13). The (64,40) subcode 
C with the best SDP is generated by G = G — {r{ 0 G1/2, r\ 0 Gi/ 2 } where r\ and r\ are the 
two rows with span [2, 15] and [3, 14] in the trellis oriented generator matrix for RM(16, 11,4). 
The SDP of the minimal trellis for C is {0,4,6,8,8,11,11,11,8,11,11,11,8,8,6,4,0} 
(<Smax,i6(C) = 11). 7 rn ax(C) = {5,6,7,9,10,11} and R(C) = {rj 0 G 1/2 , r\ 0 G 2 / 3 } where r\ 
is the row with span [5,12] in the trellis oriented generator matrix for RM(16,5,8) and 
is the row with span [4, 12] in the trellis oriented generator matrix for RM(16, 11,4). Since 
1 72(C) | = 3, one can obtain at most S parallel subtrellises in any 16-section trellis for C 
without exceeding the allowable s maXtl6 (C) of 11. Choosing C / as the subcode generated by 
G' = G — 72(C), we obtain a subcode having the SDP {0, 4, 6, 8, 6, 8, 8, 8, 5, 8, 8, 8, 6, 8, 6, 4, 0}. 
Hence we can obtain a trellis for C as the union of 8 parallel isomorphic subtrellises 
T\, 0 < i < 8. The complexity parameters of the subtrellises is presented in Table 4. 
The resulting computational complexity parameters are 

Afl(T'i) = 384, and A L (T' x ) = 22,783 (4.13) 

This gives A e //( 16) = 384 + 22, 783 + 7 = 23, 174 addition equivalent operations. 

4.4. L = 32, M = 2 and L = 64, M = 1 

When L = 32, s maXi 3 2 (C) = 12. The maximum number of parallel isomorphic sub- 
trellises possible without exceeding the allowable s maXt32 (C) = 12 in any 32-section trel- 
lis for the (64,40) subcode C is at most 4. So 4 e //( 32) > 37,476. When L — 64, 
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Smax, 64(C) = 12. Furthermore, no parallel subtrellises are possible without exceeding the 
allowable s maX)64 (C) — 12. Hence A e //(6 4) = 198,000. 


4.5. Conclusion 

By comparing the A e jf(L) obtained for L — 1,2,4,8,16,32 and 64 we find that the 
non-minimal 8-section trellis has the least effective decoding computational complexity for 
the (64,40) subcode of the Reed-Muller (64,42,8) code. This trellis has a parallel structure 
with 32 parallel isomorphic subtrellises with each subtrellis having at most 64 states. 
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Table-1: Decoding Complexity for the RM(32,16,8) Code 


L 

A(T) 

'Smax,Z/(C') 

1 

1,048,560 

1 

2 

65,599 

64 

4 

6,847 

64 

8 

4,064 

256 

16 

5,159 

256 

32 

7,995 

512 


Table 2: Complexity Parameters of 4-section trellis for (64, 38) code 


i 

0 

1 

2 

3 

4 

dim (Sm) 

0 

6 

6 

6 

0 

dim(A) 

- 

10 

10 

10 

10 

dim(Pt) 

- 

5 

5 

5 

5 

dim(<5,) 

- 

0 

6 

6 

6 

dirn(Ai) 

6 

6 

6 

0 

- 


Table 3: Complexity Parameters of 8-section trellis for (64,35) code 


i 

0 

1 

2 

3 

4 

5 

6 

7 

8 

dim(SW) 

0 

6 

6 

6 

3 

6 

6 

6 

0 

dim(.Di) 

- 

6 

6 

6 

6 

6 

6 

6 

6 

dim(-P;) 

- 

1 

1 

1 

1 

1 

1 

1 

1 

dim(£ t -) 

- 

0 

3 

3 

6 

3 

3 

3 

6 

dim(A^) 

6 

3 

3 

3 

6 

3 

3 

0 

- 


Table 4: Complexity Parameters of 16-section trellis for (64,37) code 


i 

0 

B 

B 

B 

B 

B 

B 

B 

B 

B 

10 

11 

12 



15 

16 



B 

□ 

B 

B 

B 

B 


B 

B 

8 

8 

6 

8 

6 

4 

0 

dim(D t ) 



B 

3 

3 

3 

3 

B 

B 

B 

3 

3 

3 

3 

3 

3 

3 

dim(P,-) 

B 


B 

B 

B 

B 

B 

i 

B 


0 

0 

0 

0 

0 

0 

0 

dim(^) 


B 

B 

B 

B 

B 

B 

B 

B 


3 

3 

3 

i 

3 

3 

4 

dim(A t ) 

B 

B 

B 

B 

B 


B 

B 

B 

B 



3 

1 

1 

0 

- 


15 

















































